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Abstract 



Problems of calculus of variations with variable endpoints cannot be solved without 
transversality conditions. Here, we establish such type of conditions for fractional variational 
^ ' problems with the Caputo derivative. We consider: the Bolza-type fractional variational 

problem, the fractional variational problem with a Lagrangian that may also depend on the 
unspecified end-point ip{b), where x — ip{t) is a given curve, and the infinite horizon fractional 
variational problem. 
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O I 1 Introduction 



The calculus of variations is concerned with the problem of extremizing functionals. It has many 
applications in physics, geometry, engineering, dynamics, control theory, and economics. The 
formulation of a problem of the calculus of variations requires two steps: the specification of a 
performance criterion; and then, the statement of physical constraints that should be satisfied. 
■ The basic problem is stated as follows: among all differentiable functions x : [a, 6] — K. such that 

x{a) — Xa and x(b) = Xb, with Xa, Xb fixed reals, find the ones that minimize (or maximize) the 
functional ^ 

J{x) = / L{t,x{t),x'{t))dt. 

J a 

One way to deal with this problem is to solve the second order differential equation 

dL d dL _^ 
dx dt dx' ' 

called the Euler-Lagrange equation. The two given boundary conditions provide sufficient infor- 
mation to determine the two arbitrary constants. But if there are no boundary constraints, then 
we need to impose another conditions, called the natural boundary conditions (see e.g. [TO]), 



dL 



dx' 



~ and 

t—a 



dL 



dx' 



= 0. (1) 

t=fc 



Clearly, such terminal conditions are important in models, the optimal control or decision rules 
are not unique without these conditions. 
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Fractional calculus deals with derivatives and integrals of a non-integer (real or complex) 
order. Fractional operators are non-local, therefore they are suitable for constructing models 
possessing memory effect. They found numerous applications in various fields of science and 
engineering, as diffusion process, electrical science, electrochemistry, material creep, viscoelasticity, 
mechanics, control science, electromagnetic theory, ect. Fractional calculus is now recognized as 
vital mathematical tool to model the behavior and to understand complex systems (see, e.g., [20l 
ESI |32j im |39l HH |46l [52] ) ■ Traditional Lagrangian and Hamiltonian mechanics cannot be used with 
nonconservative forces such as friction. Riewe |51j showed that fractional formalism can be used 
when treating dissipative problems. By inserting fractional derivatives into the variational integrals 
he obtained the respective fractional Euler-Lagrange equation, combining both conservative and 
nonconservative cases. Nowadays the fractional calculus of variations is a subject under strong 
research. Investigations cover problems depending on Riemann-Liouville fractional derivatives 
(see, e.g., H il [H HH HI]), the Caputo fractional derivative (see, e.g., [110111111113011111111]), 
the symmetric fractional derivative (see, e.g., [37]), the Jumarie fractional derivative (see, e.g., 
[3 [331134]), and others [H H [15l [H [Ml [13 [28] . 

The aim of this paper is to obtain transversality conditions for fractional variational problems 
with the Caputo derivative. Namely, three types of problems are considered: the first in Bolza 
form, the second with a Lagrangian depending on the unspecified end-point ip(b), where x = ip{t) 
is a given curve, and the third with infinite horizon. We note here, that from the best of our 
knowledge fractional variational problems with infinite horizon have not been considered yet, and 
this is an open research area. 

The paper is organized in the following way. Section [H presents some preliminaries needed in 
the sequel. Our main results are stated and proved in the remaining sections. In Section [3] we 
consider the Bolza-type fractional variational problem and develop the transversality conditions 
in a compact form. As corollaries, we formulate conditions appropriate to various type of variable 
terminal points. Section [4] provides the necessary optimality conditions for fractional variational 
problems with a Lagrangian that may also depend on the unspecified end-point ip{b), where x = 
ip{t) is a given curve. Finally, in Section [5] we present the transversality condition for the infinite 
horizon fractional variational problem. 

2 Preliminaries 

In this section we present a short introduction to the fractional calculus, following [26l [36l [47] . In 
the sequel, a G (0,1) and F represents the Gamma function: 

/"OO 

r(z) = / t^-^e-^dt. 
Jo 

Let / : [a, 6] — >■ M be a continuous function. Then, 

1. the left and right Riemann-Liouville fractional integrals of order a are defined by 

r(a) Ja 

and ^ 

r(a) Jx 

respectively; 

2. the left and right Riemann-Liouville fractional derivatives of order a are defined by 

aD^fix) = ^,J—^±J\x - t)-"/(t)dt. 
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and ^ 

,D^f{x)= f {t-x)-"f{t)dt, 

i (1 — aj ax 

respectively. 

Let / : [a, 6] R be a differentiable function. Then, 

1. the left and right Caputo fractional derivatives of order a are defined by 

r(i - a) J a 

and ^ 

= -FTtN / - x)-'^f'{t)dt, 
Til - a) 

respectively. 

Observe that if a goes to 1, then the operators ^D" and aD" could be replaced with ^ and the 
operators ^D^ and could be replaced with — ^ (see [17]). Moreover, we set al^f = xljjf ■= /• 
Obviously, the above defined operators are linear. If / G C^[a, 6], then the left and right Caputo 
fractional derivatives of / are continuous on [a,b] (cf. |36], Theorem 2.2). In the discussion to 
follow, we will also need the following fractional integrations by parts (see e.g. [3]): 

/' g{x) ■ ^D^f{x)dx = fix) ■ xD^g{x)dx + [jl'^gix) ■ f{x)] (2) 

J a J a 

and ^ ^ 

/ gix)-^Dmx)dx= [ fix)-aD^gix)dx-[al'y''gix)-fix)]lzl. 

J a J a 

Along the work, and following [M], we denote by diL, i = l,...,m (to g N), the partial 
derivative of function L : — )■ M with respect to its ith argument. For simplicity of notation we 
introduce operators [x] and {x, (/?} defined by 

[x]{t) = {t,x{t),^D:x{t)), 

{x,^}{t,T) ^ {t,x{t),^^D'^x{t),^{T)). 



3 Trans versality conditions I 

Let us introduce the linear space {x,t) e C^{[a,h]) x R endowed with the norm ||(a;,t)||i.oo := 
maxa<t<b \x{t)\ + maxa<f<b \'^Dfx{t)\ + \t\. 

We consider the following type of functionals: 

J{x,T)^ f L[x]{t)dt + (f){T,x{T)), (3) 

J a 

on the set 

D = {{x,t) € C\[a,b]) X [a,b]\x{a) ^ Xa} , 

where the Lagrange function L : [a, 6] x R^ — >■ R and the terminal cost function (/> : [a, 6] x R — > R 
are at least of class C^. Observe that we have a free end-point T and no constraint on x{T). 
Therefore, they become a part of the optimal choice process. We address the problem of finding 
a pair {x,T) which minimizes (or maximizes) the functional J on D, i.e., there exists (5 > such 
that J{x,T) < J{x,t) (or J{x,T) > J(x,t)) for all {x,t) e D with \\(x - x,t - T)\\i^oo < <5. 



3 



Theorem 1. Consider the functional given by ([3]). Suppose that {x,T) gives a minimum (or 
maximum) for functional ^ onD. Then x is a solution of the fractional Euler-Lagrange equation 

d2L[x]{t) + tD^{d3L[x]{t))=0 (4) 

071 the interval [a, T] and satisfies the transversality conditions 

L[x]iT)+di(biT,xiT))-x'{T) [j^-"dsL[x]it)]^^^^0 
[t4-"a3iN(i)],^y + d2^iT,x{T)) = 0. ^''^ 

Proof. Let us consider a variation {x{t) + eh{t),T + eAT), where h e C^{[a,b]), AT e M and 
e e M with |e| ^ 1. The constraint x{a) = Xa implies that all admissible variations must fulfill the 
condition h{a) = 0. Define j(-) on a neighborhood of zero by 

j{e) = J{x + eh,T + eAT) 

rT+eAT 

L[x + eh]{t) dt + 4){T + eAT, {x + eh){T + eAT)). 
If {x,T) minimizes (or maximizes) functional ([3]) on D, then j'(0) = 0. Therefore, one has 

= J [d2L[x]it)hit)+d3L[x]{t)^D^h{t)] dt + L[x]iT)AT 
A0(T, a;(T)) AT + a20(T, x{T))[h{T) + x'{T)AT]. 

Integrating by parts (cf. equation and since h{a) = 0, we get 

0=1 [a2T[x](<)+tZ??(a3iN(t))]/i(t)rfi+ [t4-"(a3T[x](<))/i(t)]^^^ + L[x](T)AT 

+\(biT, x{T))AT + d2(l){T, x{T))[h{T) + x'{T)AT] 

[d2L[x\{t) + tD^{d:iL[x\{t))] hit) dt 

+AT [L[x\{T) + di^{T,x{T))~x'{T) [tl]r--d3L[x\{t)]^^^] 
+ [[t4-"a3L[a;](t)]^^j, + 920(r,x(T))] [h{T) + x' [T) AT] . 

As h and AT are arbitrary we can choose h{T) = and AT = 0. Then, by the fundamental 
lemma of the calculus of variations we deduce equation (j4]) . But if x is a solution of (jl]) , then the 
condition ([6|) takes the form 

= AT [L[x\{T) + d,4>{T,x{T)) - x'{T) [a]r<^d:iL[x]{t)]^^^] 
+ [[t4-"93iNW]t^T + ^20(r,x(T))] [h{T)+x'{T)AT]. 

Restricting ourselves to those h for which h{T) — —x'{T)AT we get the first equation of ([S]). 
Analogously, considering those variations for which AT = we get the second equation of ([5]). □ 

Example 1. Let 

J{x, T)= [t^^l + i^D^xit))^] dt, T e [0, 10]. 
Jo 

It easy to verify that a constant function x{t) — K and the end-point T — 1 satisfies the necessary 
conditions of optimality of Theorem [7J with the value of K being determined by the initial-point 
x(0). 

In the case when a goes to 1, by Theorem [1] we obtain the following result. 
Corollary 1. (!22j. Theorem 2.24) ^/(2;,T) gives a minimum (or maximum) for 

J{x,T)^ [ L{t,x{t),x'{t))dt + (l){T,x{T)) 
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on the set 

{x e C^{[a,b]) I x{a) = Xa} x [a,b], 
then X is a solution of the Euler-Lagrange equation 

d2L{t,x{t),x{t)) - ^dzL{t,x{t),x{t)) = 
at 

on the interval [a,T] and satisfies the transversality conditions 

L{T, x{T),x'{T)) + 9i0(r, x(T)) - x'{T)d^L{T, x{T),x'{T)) = 0, 
d^LiT, x(T), x'(T)) + d2<j){T, x{T)) = 0. 

Now we shall rewrite the necessary conditions ([S]) in terms of the increment on time AT and 
on the consequent increment on x, Axt- Let us fix e = 1 and consider variation functions h 
satisfying the additional condition h'{T) = 0. Define the total increment by 

AxT = (x + h){T + AT) - x{T). 

Doing Taylor's expansion up to first order, for a small AT, we have 

{x + h){T+ AT) -{x + h){T) = x'{T)AT + 0{ATf , 

and so we can write h{T) in terms of AT and Axt'- 

h{T) = AxT-x'{T)AT + 0{ATf. (8) 

If {x, T) gives an extremum for J, then 

d2L[x\{t) + tD^{d:iL[x]{t)) = holds for all t e [a, T]. 

Therefore, substituting ([8]) into equation (|6]) we obtain 

AT [L[x]{T) + ai0(T, x{T)) - x'{T) [j^^-^dMit)],^^] 

+ AxT[[t4~"53iNW]t=T + ^20(r,x(T))] +O(AT)2 = 0. (9) 

We remark that the above equation is evaluated in one single point x ^ T. Equation ^ replaces 
the missing terminal condition in the problem. 

Let us consider five particular cases of constraints that can be specified in the optimization 
problem. 

A. Vertical terminal line 

In this case the upper bound T is fixed and consequently the variation AT = 0. Therefore, 
equation ([9]) becomes 

AxT [[tI^-''d3L[x]{t)]^^^ + d2HT,x{T))] = 0, 
and by the arbitrariness of Axt, we deduce 

[tI^-''d3L[x]{t)]^^^ + d2^{T,x{T)) = 0. 
When = 0, we get the natural boundary condition as obtained in [T], equation (25): 

[tl].-'^d^L[x]{t)\^^^ = Q. 

For = and a H> 1, we get the second equation of ([T]). 

B. Horizontal terminal line 
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(10) 



Now we have Axt — but AT is arbitrary. Hence, equation implies 

L[x]{T) + di^{T,x{T)) - x'{T) [j^-'^dsLixKt)]^^^ = 0. 

When = and a — > 1, we get equation (3.11) of It is worth pointing out that the integer 
case has an economic interpretation (see explanation in [23] . pags. 63-64). 

C. Terminal curve 

In this case the terminal point is described by a given curve ip, in the sense that x(T) = ijjiT), 
where ip : [a. 6] ^ R is a prescribed differentiable curve. For a small AT, from Taylor's formula, 
one has 

AxT = ^{T + AT) - ijj{T) 
= 'ijyiT)AT + 0{AT) 

Substituting ^ into dH) yields 

(V'(r) - x'{T)) { [tI^~"dsL[x]{t)]^^^ + 920(T,x(T))} + L[x]{T) + di<l>{T,x{T)) = 0. 

For = 0, we obtain equation (29) of [2 . For (p = and a — > 1, we have equation (3.12) of |23) . 

D. Truncated vertical terminal line 

Now we consider the case where AT = and x(T) > Xmin- Here Xmin is a minimum permissible 
level of X. By the Kuhn- Tucker conditions, we obtain 

[tI^^''d3L[x]{t)]^^^ + d2(b{T,x{T)) < 0, X{T) > Xrmn, 

{x{T) - x™„) { [t4-"a3iN(t)]^^^ + d2^{T,x{T))) = 
for maximization problem; and 

[t4-"a3iN(i)]t^T. + d2^{T,x{T)) > 0, x{T) > x™„, 

{x{T) - x™„) { [tI^-'^d3L[x]it)]^^^ + a2</)(T,a:(T))} = 

for minimization problem, li = and a 1, we get equations (3.17) and (3.17') of j23j. 

E. Truncated horizontal terminal line 

In this situation, the constraints are Axt — and T < T^ax- Therefore, as in the previous 
case, we obtain: 

L[x\{T)+di4>{T,x{T)) - x'{T) [tl]r''d3L[x]{t)\^^^ > 0, T< Tma., 

{T - T„,ax) [L[x]{T) + di<P{T,x{T)) - x' {T) [tl]r'-d3L[x]{t)\^^^] = 
for maximization problem; and 

L[x\{T)+di4>{T,x{T)) - x'(T) [tI^-''d3L[x]{t)]^^^ < 0, T< T^ax, 

{T-T,na.) [L[x]{T)+dMT,x{T))-x'{T) [tI^-''d3L[x]{t)]^^^] =0 

for minimization problem. Observe that for (/) = and a — > 1, we get equation (3.18) and (3.18') 
of [23]. 

4 Transversality conditions II 

In this section we consider the following variational problem: 

Jix,T) 



I L{x, ip)(t, T)dt — > extr 

. ^ (11) 

x{T) = cp{T) 



where L : [a, 6] x R'^ -> K and (p : [a, fe] — > R are at least of class C^. Here x — (p{t) is a specified 
curve. 
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Theorem 2. Suppose that {x,T) is a solution to problem (jlip . Then x is a solution of the 
fractional Euler- Lagrange equation 

d2L{x, ifijit, T) + tD^id^x, T)) = (12) 

071 the interval [a,T], and satisfies the transversality condition 

{^'{T)-x'{T))[tI^-"d3L{x,^}{t,T)]t^T+^'{T) f diL{x,^}it,T)dt+L{x,^}{T,T) = 0. (13) 

Proof. Suppose that (x, T) is a solution to problem (jlip and consider the value of J at an admissible 
variation {x{t) + eh{t), T + eAT), where e e M is a small parameter, AT e M, and h E C^{[a, b]) 
with h{a) = 0. Let 

pT+eAT 

j{e)^J{x + eh,T + eAT)= / L{x + €h,ip}{t,T + eAT)dt. 

J a 

Then, a necessary condition for {x, T) to be a solution to problem (|lip is given by 

j'(0) = ^ j [d2L{x, ^}(t, T)h{t) + a3L{x, T)'^D^h(t) + d^Hx, <f}{t, T)^'{T)AT] dt 
+L{x,<f}{T,T)AT^O. 

Integrating by parts (cf. equation ©), and since ft, (a) = 0, we get 



[d2L{x,^}{t,T) + tD^{d3L{x,^}it,T))]h{t)dt+[tI^-''{d3L{x,ip}{t,T))h{t)] 

+L{x, ip}{T, T)AT + /J [diL{x, if}{t, T)ip'{T)AT] dt. 



t=T 



(14) 

As h and AT are arbitrary, first we consider h and AT such that h{T) = and AT = 0. Then, by 
the fundamental lemma of the calculus of variations we deduce equation (jl2|) . Therefore, in order 
for (x, T) to be a solution to problem (fTTj) . x must be a solution of the fractional Euler-Lagrange 
equation. But if x is a solution of (fT2|) . the first integral in expression (fT4| vanishes, and then the 
condition ([T4|) takes the form 

[tI^-"{d3L{x,if}it,T))hit)]^^^+L{x,ip}iT,T)AT+ / [diL{x,ip}it,T)ip'{T)AT]dt ^ 0. (15) 

^ a 

Since the right hand point of x lies on the curve z — ip{t), we have 

x{T + eAT) + eh{T + eAT) = ip{T + eAT). 
Hence, differentiating with respect e and setting e = we get 

hiT) = {ip'{T) ~ x'{T))AT. (16) 

Substituting HH) into ^ yields 

AT{ip'{T)^x'{T)) [j]r^d:,L{x,ip}{t,T)]^^^ 

+ ip\T)AT I d4L{x, ip}{t, T)dt + L{x, ip}{T, T)AT = 0. 

Since AT can take any value, we obtain condition (|13p . □ 
In the case when a goes to 1, by Theorem [2] we obtain the following result. 
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Corollary 2. If {x,T) gives an extremum for 

J{x,T)^f L{t,x{t),x'{t),if{T))dt 

J a 

on the set 

{x e C^{[a,b]) \x{a) = Xa} x [a,b], 
then X is a solution of the Euler-Lagrange equation 

d2Lit,xit),x'it),ip{T))- ■^^d3L{t,x{t),x'it),ip{T)) ^0, for all te [a,T], 
and satisfies the following transversality condition 

(^'(T) - x'{T))dsL{T, x{T),x'{T),ip{T)) + ^\T) f diL{t, x{t),x'{t), ^(T)) dt 

+L{T,x{T),x'{T),cpiT))^0. 
Example 2. Consider the following problem 

J(x, T)= [ [{^D?xit)y + ^2(T)] dt — > min 

(17) 

x{0) = 1. 

For this problem, the fractional Euler-Lagrange equation and the transversality condition (see 
TheoremW^ are given, respectively, by 

tD^ i^D-^xit)) = 0, (18) 

2(1 - x'{T)) [j^-'^^D'^xit)] + + {^D'^x{T)f +T' = 0. (19) 

Note that it is a difficult task to solve the above fractional equations. For < a < 1 a numerical 
or direct method should be used jlOf - When a goes to 1, problem (jl7p becomes 



J(x,T)= r [x\t) f + ^\T)\dt 
Jo 

x{T) = ^(T) = T, 
a;(0) = 1. 



(20) 



and equations (|18p - (|19p are replaced by 

x"{t) - 0, (21) 

2(1 - x'{T))x'{T) + 2T^ + (.t'(T))2 + = 0. (22) 
Solving equations ()2ip and ()22p we obtain that 



i(0 = ^^^^t+l, r= 1^-6 + 6713 



is a candidate for minimizer to problem ((20)) . 
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5 Infinite horizon fractional variational problems 



Starting with the Ramsey pioneering work |50) . infinite horizon variational problems have been 
widely used in economics (see, e.g., [T51 [531 [23 HOI Ei] and the references therein). One may assume 
that, due to some constraints of economical nature, the infinite horizon variational problem does 
not depend on the usual derivative but on the left Caputo fractional derivative. In this condition 
one has to consider the following problem: 



r+oo 

J(x) = / L[x]{t)dt 



(23) 



where L : [a, +oo] x ^ R is at least of class and 

Doo = {x e C^{[a, +oo]) I x{a) = Xa} ■ 

The integral in problem (1231) may not converge. In the case where the integral diverges, there may 
exist more than one path that yields an infinite value for the objective functional and it would be 
difficult to determine which among these paths is optimal. To handle this and similar situations 
in a rigorous way, several alternative definitions of optimality for problems with unbounded time 
domain have been proposed in the literature (see, e.g., [HI [33 13S] ) • Here, we follow Brock's notion 
of optimality [TH], i.e., a function x G Doo is said to be weakly maximal to problem (153)) if 



T 



lim inf / \L\x](t) - L\x](t)] dt<0, 

-*+ooT'>T 



for aU X e Doo ■ 

Using similar approach as in |40i 145] , we obtain the transversality condition for infinite horizon 
fractional variational problems. First, a lemma that we will use in the proof of Theorem [31 



Lemma 1. If g is continuous on [a, +oo) and 

r' 

im inf / 

T- 

for all continuous functions h : [a, +oo) — > R with h(a) = 0, then g{t) — 0, for all t > a. 



lim inf / g(t)h(t) dt = 

^^+oo T'>T V ; 



Proof. Can be done in a similar way as the proof of the standard fundamental lemma of the 
calculus of variations (see, e.g., [19 J. □ 

Theorem 3. Suppose that x is a weakly maximal to problem (|23p. Let h e C^([a, +cx)]) and 
e e M. Define 

Aie,T') . r L[x + eh]it)-L[x]it)^^,^ 
Ja e 

f-T' 



V{e,T) = ^nf^^ [L[x + eh]{t) - L[x]{t)]dt; 



V{e) = hm V{e,T). 



Suppose that 



1. lim — ^— — ^ exists for all T; 
e^O e 

2. lim — — — ^ exists uniformly for e; 
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3. For every T' > a, T > a, and e G R \ {0}, there exists a sequence {A{e,T^))^^j^ such that 

lim A(e,T^)= inf A(e,T') 

n— >+oo T'>T 

uniformly for e. 

Then x is a solution of the fractional Euler- Lagrange equation 

d2L[x\(t) + tD'^{d:iL[x]{t))^G 
for all t G [a, +cxd), and for all T > t. Moreover it satisfies the transversality condition 

lim inf t47"53i[x]m]t=T'a;(T') = 0. 

Proof Observe that V(e) < for every e e M, and V{0) = 0. Thus V'(0) = 0, i.e., 

= limM 



Km inf / \d2L\xMt)h(t)+d-iL\x]{t)^D?h{t)\dt. 

T^+ooT'>T '- I Jwa I V /J 

Thus, integrating by parts and since h{a) = 0, we get 

^lim^^nf^ 1^ [d2L[x\{t) + tD'^,{d:iL[x]{t))]h{t)dt + [tI^^7'^d:iL[x\{t)h{t)]t=T' | = 0. (24) 
By the arbitrariness of ft, we may assume that h{T') — and so 

r 

T^+oo T'>T 

Applying Lemma [1] we obtain 

d2L[x\{t) + tD^{d^L[x\{t))^Q 
for all t e [a, +oo), and for all T > t. Substituting equation (|25l) into equation ([M)) we get 

lim inf [,47"93iN(i)]t=T'/i(T') = 0. (26) 

To eliminate /i from equation (1261) . consider the particular case h{t) — x(i)a;(i), where x '■ 
[a, +oo) R is a function of class such that x(a) = and x{t) = const ^ 0, for all t > to, for 
some to > a. We deduce that 



\hn^M^£ [d2L[x]{t) + tD^,id3L[x]mh{t)dt - 0. (25) 



lim inf [tI^7"d3L[x]{t)]t=T'x{T') = 0. 

T->+ooT'>T ^ 



□ 



6 Conclusion 

Transversality conditions are optimality conditions that are used along with Euler-Lagrange equa- 
tions in order to find the optimal paths (plans, programs, trajectories, etc) of dynamical models. 
The importance of such conditions is well known in economics models and other phenomena 
whose effects can be spread along time, e.g., radioactive, pollution. In this paper we have given 
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in the compact form transversality conditions for fractional variational problems with the Caputo 
derivative. The fractional variational theory is only 15 years old so there are many unanswered 
questions. For instance, the question of existence of solutions to fractional variational problems 
is a complete open area of research (more about this important issue the reader can find in 
[m |43l [53l [54] ) . Other interesting open question is about the convergence of the objective func- 
tional J{x) = 1'!'°° L[x]{t)dt. This issue should be treat carefully, especially as the functional 
depends on fractional type of derivatives. Generally, in order to solve fractional Euler-Lagrange 
differential equations and apply transversality conditions one needs to use numerical methods. 
Some progress is being made already on the numerical methods for fractional variational problems 
m m im [49] but there is still a long way to go. We believe that aU pointed issues will be considered 
in a forthcoming papers. 
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